Motion of nanodroplets near chemical heterogeneities 
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We investigate the dynamics of nanoscale droplets in the vicinity of chemical steps which separate 
parts of a substrate with different wettabilities. Due to long-ranged dispersion forces, nanodroplets 
positioned on one side of the step perceive the different character of the other side even at some 
distances from the step, leading to a dynamic response. The direction of the ensuing motion of such 
droplets does not only depend on the difference between the equilibrium contact angles on these two 
parts but in particular on the difference between the corresponding Hamaker constants. Therefore 
the motion is not necessarily directed towards the more wettable side and can also be different from 
that of droplets which span the step. 



The interest in dynamical wetting phenomena has in- 
creased significantly with the development of micro- and 
nanofluidic systems, which allow one to handle minute 
amounts of liquid containing, e.g., DNA or proteins for 
chemical analysis and biotechnology '5] . In particular 
chemically heterogeneous systems with tailored, spatially 
varying wetting properties have found important applica- 
tions in this context [3, 0|- open microfluidic systems, 
fluids are transported on chemical channels, i.e., lyophilic 
stripes embedded in lyophobic substrates. While present 
devices are based mostly on micron sized channels, fur- 
ther miniaturization down to the nanoscale is clearly on 
the roadmap. This will eventually lead to nanofluidic sys- 
tems for which a variety of physical phenomena, which on 
the microscale and above are either irrelevant or summa- 
rized into boimdary conditions, become important 0, Q. 
For the optimization of the performance of nanofluidic 
systems it is critical to understand the basic fluidic issues 
occurring on those scales. Recent theoretical studies of 
nanoscale fluids on chemically 7] and topographically Q 
structured substrates have underscored the importance 
of such investigations. 

Those analytical tools [1, [l^, [ll|, which rely on clas- 
sical macroscopic theory are not adequate for this pur- 
pose. Apart from some molecular dynamics simulations 
01 13 1 which are computationally demanding, to a large 
extent the available numerical inve stig ations are based 
on solving thin fllm equations [ij, [3 dj [13, EBl- In 
most of theses studies the chemical heterogeneities are 
introduced via abrupt, lateral changes of those param- 
eters which characterize the potentials of homogenous 
substrates [13, [H, [l3|- However, this does not capture 
the actual behavior of such substrate potentials, even if 
the underlying chemical steps are taken to be atomically 
sharp [l^. Smooth chemical heterogeneities have been 
studied in Ref. [l8| by introducing a continuously vary- 
ing Hamaker constant. But this approach is only applica- 
ble for very smooth variations of the wetting properties. 
Studies of the dynamics of droplets in the vicinity of to- 
pographic steps have shown that on the nanoscale a de- 
tailed modeling of the substrate and thus of the resulting 
effective interface potential is mandatory . For a chem- 



ical step, this has been worked out within the framework 
of microscopic density functional theory with a view on 
the morphology of static wetting films [l^, [U [23] . 

Here we study the dynamics of nanodroplets in the 
vicinity of chemical steps, i.e., sharp and straight bound- 
aries between two parts of a substrate with different 
wetting properties. Such steps appear in open micro- 
and nanofluidic systems as the edges of chemical chan- 
nels. As a paradigmatic case we focus on a chemical 
step formed by two adjacent quarter spaces composed 
of different substrate particles. We analyze the driving 
force on the droplets and perform numerical calculations 




FIG. 1: The chemical step is translationally invariant along 
the z axis (orthogonal to the image plane). A nanodroplet 
(full line) is exposed to the vertically and laterally varying 
DJP, the contour plot of which is shown. Here both sides 
of the substrate are taken to correspond to the Q case with 
= 6^ and = 97.2° {B^ = 0, = 3) and 9% = 51.3° 
(B^ = 0, = 1) [see Eqs. Q and © and the main text]. 
Lengths are measured in units of 6^ (see the main text). A 
coating of the substrate is not indicated as here B"^ = B^ = 
so that for x ±oo the equilibrium thicknesses of the 
underlying film are the same, i.e., — Vq = yo = . 



2 



employing the standard boundary integral method for 
hydrodynamic Stokes flow providing the underlying dy- 
namics. 

As depicted in Fig. [1] we consider a partially wet- 
ting, non-volatile, and incompressible liquid forming 
a nanodroplet on top of a precursor wetting film over 
a chemical step. The type of chemical step which 
we consider here can be viewed as being composed 
of two different quarter spaces with each of the two 
corresponding upper surfaces coated additionally with 
a thin layer of different materials. Accordingly, as a 
basic element we first consider an edge without surface 
coating, say the left part of the substrate shown in 
Fig. [T] which we denote by the superscript <] . Assuming 
Lennard-Jones type intermolecular pair potentials 
Vafsir) = Mafs/r^"^ — Nap/r^, where Mafs and Naf3 are 
material parameters, and a and /? relate to liquid (I), 
substrate (s), or coating (c) particles, the local disjoining 
pressure (DJP) corresponding to an edge occupying 



tension. Thus the dimensionless DJP 11*^ 



nf = {r e : 



n:(x,y) 



< 0, y < 0, z e M} is given by |2S 



n< r 



r - 



d'r', (1) 



with AM<= = pfMu - pipfMfi and AiV<= = pjNu - 
pipfN^i where pi and pf are the number densities of the 
liquid and the substrate, respectively. Because of its low 
density in Eq. ([T|) the effect of the vapor or gas phase 
has been neglected. The edge geometry allows one to 
analytically calculate both integrals in Eq. Ilj which we 
denote as Ilg^^ and Ilg^, respectively, so that 11^ (x, y) = 

Likewise, the contribution to the DJP of a thin coating 
layer of thickness on the upper part of an edge 
f^^^ = {r e R3 I 2; < < y < 0,2 G M} can 

be determined by assuming also a van der Waals type 
interaction between the coating and the liquid parti- 
cles, i.e., n^^(a;,?/) = - AiV'^/|r - r'|^ d^^', with 
AN'^ = pfNu - pip^Nf and U^<{x -00, y) = 
—TrAN'^d'^/{2y'^) As a simplification we have ne- 
glected the effect of the repulsive part of the liquid- 
coating interaction which gives rise to a contribution 

shorter ranged (~ y"^*^) than the corresponding one 

ni2<i ^ y-9 [23, Hj. The DJP of an edge including 
the coating of its upper side is U^^{x,y) = IIf{x,y + 
d'^) + ir^'^{x,y). For x —00 and much smaller 
than the thickness of the wetting layer y^, the DJP 
of the coated edge reduces to that of a coated, later- 
ally /lomogeneous substrate: n5j(y) — ttAA/ ^/(45y^) — 
7rAiV^/(6y3) + 7rA7V"^d^/(2y''), with AA^"^ = AA^^ - 
AN'< = pi{p<N< - pfN<); U<^{y = y^) = 0. 

At this point we introduce dimensionless quantities 
(marked by a star) such that lengths are measured 
in units of &^ = [2AA^^/(15|AAf<=|)]^/^ which for 
AA^^ > is the equilibrium wetting film thickness on 
the uncoated homogeneous substrate. The DJP is mea- 
sured in units of where a is the liquid- vapor surface 



from the edge {x — > —00) has the following form: 



(2) 



in the following we drop the stars. 

In Eq. (0) = TrAN"^d^/{2A-^b^'^) quantifies the 
contribution of the coating and — A'^h'^/a where 

A-^ = 7r(AM^/45)"^/^(|A7V^|/6)^/^ measures the ef- 
fect of the competing intermolecular forces relative to 
the surface tension of the liquid-vapor interface. Since a 
more refined analysis of the DJP beyond Eg. ([T]) yields 
B ^ even in the absence of a coating layer [2J, [25| , in 
the following we consider B as an independent disposable 
parameter. In the second term on the right hand side of 
Eq. ^ the upper plus (lower minus) sign corresponds to 
AA^<' < (AN^ > 0). In the following we shaU refer to 
these cases as the plus © and the minus G cases, respec- 
tively. The dimensionless form of the DJP (in units of 
cr/&^) for a single edge coated on the upper side is given 

by 



Il<{x,y) ^ 



45Ill^<ix,y) ^6nl<{x,y) 



TT AM<i 
2B< [^n^<{x,y)] 
ttAN'o 



± 



|AA^<| 



(3) 



Here 11^^^, Ilg^, and 11"^ are measured in units of a/b^ 
whereas AM^, AA^^, and AA^'^ are taken in units of 
(t(6^)*, cr(6^)^, and (j{b^)^, respectively; x and y are 
measured in units of 6^ so that and B"^ are dimen- 
sionless. The contact angle 6^^ of macroscopic droplets 
and those values of and B'^, which give rise to par- 
tial wetting, are related via cos0^ = 1 + Il'^f^{y) dy 



with y^ given imphcitly by U-^i^iy^) = [8|, |25| so that 
yQ{B^ = 0) = 6^. Within these admissible ranges of val- 
ues for and B^ (see the insets of Fig. 2 in Ref. [8]), 
the actual contact angle of nanodroplets defined, e.g., 
via their slope of the point of inflection or by spheri- 
cal extrapolation of their top cap towards the substrate, 
may differ from that of macroscopic droplets depending 
on the size of the nanodroplets and details of the DJP. 
With q = b^ /b"^ the dimensionless form of the DJP for 
the right edge is given by 



^, A5q^Ul^<{-x,y) 



± 



6q^Ul<i-x,y) 



TT AM< 

2q^B> [-U"^<i~x,y)] 
ttAN'< 



TrlAA^^ai 



(4) 



Due to the additivity of the interatomic potentials used 
here, the DJP of the chemical step can be obtained by 
superimposing the DJP of the constitutional parts, i.e., 
the two edges coated on the upper side. Thus the DJP 
of the whole substrate with a chemical step is given by 
U{x,y) = Il^^{x,y) + U^^{x,y). Figure [1] shows the con- 
tour lines of the DJP across a chemical step for a typical 
example. 
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FIG. 2; The motion of a nanodroplet across three different 
chemical steps for the case on both sides and for = 0, 
B'^ = 0, and q = /h^ = 1.0 so that = = j/o- 
The uppermost part shows the initial interface profile. The 
parameters for the lower three parts are [ 6*^^ = 97.2° (C^ = 
3), e% = 51.3°(C^ = 1)], [6^^ = 97.2°(C<^ = 3), 0% = 
75.5° (C^ = 2)], and [6*,^, = 120° (C^ = 4), 6% = 75.5° (C^ = 
2)] from top to bottom. The profiles correspond to times 
t = 355, 2000, and 22025; t = 375, 8125, and 24300; and 
t = 185, 7800, and 23500; respectively, in units of ^6^/(C^cr). 
< (>) is the less (more) wettable substrate. 



In order to probe the influence of the DJP on droplets 
near and on chemical steps we have performed meso- 
scopic hydrodynamic calculations based on the two- 
dimensional Stokes equation. In dimensionless form 
the continuity and Stokes equation read V ■ u = 
and V^u = Vp, respectively, where u = {ux, Uy) 
is the velocity vector and p is the hydrostatic pres- 
sure. With the viscosity /i, the velocity and time scales 
are taken to be C^cr//i and fib'^ /{C^a), respectively. 
Lengths and pressure are expressed in units of 6^ and 
cr/&^, respectively. In the above Stokes equation the 
factor appears because its dimensionless form has 
been obtained by rescaling with the parameters of the left 
hand side of the substrate. We have solved theses equa- 
tions numerically with a standard biharmonic boundary 
integral method [2^. A no-slip boundary condition has 
been employed for the impermeable liquid-solid interface 
and it has been imposed that there is no flux through 
the end sides of the system. Along the liquid-vapor in- 
terface it has been assumed that the tangential stresses 
are zero and that normal stresses are balanced by pres- 
sure, surface tension, and the disjoining pressure, i.e., 
n • r • n = —p + an + 11 with the local curvature k, the 
stress tensor r, and the unit surface normal vector n. 
The initial droplet shape has been taken to be a hemi- 
sphere which is smoothly connected to the wetting layer, 
i.e., y{x;t = 0) = y< + a{l - - g)/a]^]}\\-\-3\-+i 
with the droplet height a in the center equal to half the 
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FIG. 3: (a) 9{x) and (b) dimensionless lateral force densities 
f{x) for droplets with a = 15 (for x see Fig. [5)) in the vicinity 
of a chemical step separating two parts of the substrate with 
= 3, = = 0, and (7=1. Both substrate sides 
correspond to the Q case. / is measured in units of cr/(6^) . 



base width and the distance g of the droplet center from 
the boundary (see Fig. [5]). In this study m was chosen 
to be 10. 

We first consider the case without coating (B^ = 
= 0, which implies that AA^ is negative on both 
sides corresponding to the case). Numerical results for 
different values of and are shown in Fig. [21 At 
time t = droplets of height a = 15 have been posi- 
tioned at a distance w = 5 (see Fig. [2]) from the chemical 
step on the less wettable substrate. was selected such 
that, during the initial relaxation of the prepared droplet 
shape, w does not decrease. Thereafter, in all cases shown 
in Fig. [2] the droplets gradually move towards the more 
wettable side (>), as expected intuitively, and continue 
their motion there. This indicates that the nanodroplets 
can perceive the presence of the other part of the sub- 
strate over some lateral distances. As a function of time 
the wetting layer thickness changes slightly which is ex- 
pected due to the Laplace pressure in the droplet and in 
view of the lateral boundary conditions on the flux. 

On a homogeneous substrate, the equilibrium con- 
tact angle can be calculated from the DJP. Extending 
this formula, one can define a spatially varying "con- 
tact angle" via cos6'(a:;) = 1-1- J^i^^^U{x,y)dy, with 
n(x,yo(a;)) = 0. For n(x,y) rapidly varying as a 
function of a; (e.g., in the close vicinity of a chemical 
step) and for very small droplets, 9{x) is not the ac- 
tual contact angle, but nonetheless, as intuitively ex- 
pected, dxO{x) is related to the lateral DJP induced 
force density (in units of (t/(6^)^) acting on a droplet: 
/ = l/r^d X ^Q^^Yi{x,y)nxdS where dVld and Vld are 
the dimensionless droplet surface and volume, respec- 
tively, and Ux is the x-component of the outward surface 
normal n. For the present purposes and for the sake of 
simplicity we have focused on nanodroplets large enough 
such that j.^^ n(a;, y)c??/ w J^f^_^^Il{x,y)dy. This means 
that for the types of the DJP considered in this study 
a is taken to be at least 10 times the thickness of the 
underlying wetting film. Within this range, the phe- 
nomena discussed below are largely independent of the 
droplet size. Since in the cases studied here the differ- 
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FIG. 4: (a) 9{x) and (b) dimensionless lateral force densities 
f{x) for droplets with a = 15 in the vicinity of a chemical 
step between two substrates with = 3, — = 0, 
and q = 1.25. Both parts of the substrate correspond to the 
case. / is measured in units of cj/{h^) . 



ences in wettability between the two parts of the sub- 
strate are relatively small, in the following we have esti- 
mated /(x) by considering droplets with the shape used 
as initial condition for the Stokes dynamics but centered 
aX X = X. For large and symmetric droplets and if 
n(a;, y) varies slowly over the width of the droplet such 
that n(a;o + a,y) - Tl{xo - a,y) w 2adxTl{xQ,y), f and 
—dO/dx have the same sign. Thus droplets will move 
towards smaller 0{x), i.e., more wettable regions. 

Figure [31^ a) shows 9(x) in the vicinity of the chemical 
step for the parameters used in Fig. O 0{x) monotoni- 
cally decreases from 6*^ to 0^ and, as shown in Fig.[3jb), 
the force acting on droplets of initial height a = 15 is 
positive for all positions x of the center of the drop. The 
force curve has a rather flat, plateau- like shape with a 
maximum at x = and varies sharply if one of the con- 
tact lines reaches the chemical step. This result is in 
agreement with the numerical calculations and indicates 
that the droplets never come to a complete stop. 

In the previous example the local contact angle 0{x) 
changes monotonically from 9^^ to 9^^. This is not neces- 
sarily the case, in particular not for steps between mate- 
rials generating different thicknesses of the wetting layer, 
and for two substrates with different forms of the DJP, 
e.g., on the left and © on the right side. 

Figure [4] shows 9{x) and f{x) for = = and 
q — 1.25. For < 1.5, 9{x) decreases monotonically 
and f{x) is positive for all x. Therefore a droplet will 
move towards the more wettable side. For = 1.6 
and larger, 9{x) is nonmonotonic. The cases between 
C'^ = 1.6 and 2.15 are particularly interesting because, 
although 9^^ > 9^^, dx9 > outside a small region 
around the step. A droplet is therefore expected to move 
towards the less wettable substrate. This is indeed the 
case because / also changes sign. Zeroes of f{x) with 
dxf{x) < indicate that the chemical step can act as a 
pinning site for droplets. For C"^ > 2.4 the two parts 
of the substrate exchange roles in that 0^ < 9^^. Al- 
most everywhere f{x) < except for those values of x 
for which one of the contact lines touches the step. Thus, 
a droplet is expected to move towards the more wettable 
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FIG. 5: Motion of a nanodroplet near a chemical step in the 
e case on both sides, ^ ^ 0, = 3, = 2.15, 
and <j = 1.25 (implying 9^^ = 97.2° and 61^, = 90.4°; compare 
Fig-H]), starting at different positions. The droplet shapes are 
shown, top to bottom, at t = 170, 100, and 45 (dashed) and 
at t = 4750, 4300, and 40000 (soUd), respectively. 



side. However, we expect pinning of the contact lines at 
the step as the droplet moves from [> to O because in 
this case dxf{x) < at the zeroes of f{x). 

In order to confirm these theoretical predictions we 
have performed a series of numerical calculations by po- 
sitioning a nanodroplet at a distance w = 5 on the right 
hand side of the chemical step. For the case = 2.15, 
i.e., if the nanodroplet is initially positioned on the more 
wettable part of the substrate, the results of these calcu- 
lations are shown in Fig. [5] for a nanodroplet with a = 15. 
This droplet moves to the left (i.e., counterintuitively to- 
wards the less wettable substrate), slows down, and stops 
when its advancing contact line touches the chemical 
step. A droplet initially spanning the step moves to the 
right, as expected, but stops when its receding contact 
line touching the step. Droplets completely located on 
the left part of the substrate move away from the chem- 
ical step, i.e., away from the more wettable substrate, 
with the velocity decreasing as the distance from the step 
increases. However, we do not observe a complete stop. 
These numerical results are in complete agreement with 
our analysis of the effective DJP induced force. Figure [H] 
shows an example for q = 1 in which the droplets can 
move towards the less wettable substrate. This is the case 
for different signs of the long-ranged term (here © case on 
the left and © case on the right) and for non-zero values 
of B (here = -2.762, = 2, and B'> = -1). These 
values are chosen such that the equilibrium wetting film 
thicknesses are equal on both sides of the substrate. Out- 
side a narrow region around the chemical step dx9{x) is 
positive everywhere. This indicates that the droplet will 
move to the left independent of the relative values of 0^ 
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FIG. 6: (a) 9{x) and (b) dimensionless lateral force densities 
f{x) for droplets with a = 15 in the vicinity of a chemical 
step separating two substrates (® case on the left and case 
on the right) with B'^ = -2.762, = 2, = -1, and 
q = 1 for different values of . f is measured in units 
of a/{b'^) . The different coatings of the two parts of the 
substrate associated with the non-zero values of and B^ 
are indicated in the figures. 

and 9^^, although it will be pinned at the step if coming 
from the right. This is confirmed by the analysis of f{x) 
in Fig. [H^b). If both contact lines are on the same side 
of the step, > a) < 0. If the droplets span the step, 
i.e., < a, the direction of the motion is towards the 
more wettable part of the substrate. 

In all the cases we considered the direction of motion 
of a droplet near a chemical step (i.e., with both con- 
tact lines on the same side of the step) is determined by 
dx0{x). The sign of 8^9 {x) is the same as the sign of 
—dxTl{x,yQ{x)). Expanding Il{x,y) for large \x\ up to 
0(|a;|^'^) one finds on both sides of the step the limiting 
value for a homogeneous substrate. This means that the 
wetting film thickness yo is independent of x up to this 
order. To leading order the gradient of the equilibrium 
contact angle with respect to x is the same on both sides 
of the chemical step and given by 

d.e{x) ^ [-{±q') C> ± C^] + 0{\x\-'), (5) 

with the plus and minus signs in front of C'^ and 
corresponding to the © and the case on the right and 
the left hand side, respectively. This means that asymp- 
totically far from the step the droplets move in the same 
direction on both sides of the step. For instance, for the 
case discussed in Fig. |4] (with B'^ = = 0) the equi- 
librium contact angle can be shown to be larger on the 



right side of the step if gC"> > C^, i.e., for > 2.4 and 
the droplets move to the left for C"> > C^/q^ = 1.536 
and \x\ > a, in agreement with the observations. 

In summary, we have outlined an approach which al- 
lows one to study in detail the behavior of nanodroplets 
near as well as on chemical heterogeneties. Our in- 
vestigation reveals the dynamics of nanodroplets in the 
vicinity of chemical heterogenities caused by long-ranged 
forces. We have shown that the direction of motion of 
the droplets is, to leading order in the distance from 
the step, determined by the competition of the van der 
Waals forces acting between the droplet and the two dif- 
ferent materials of the substrate, i.e., the difference in the 
Hamaker constants (see Eq. ([5])), rather than the equi- 
librium contact angles which depend also on the short- 
ranged parts of interaction potentials and on the sub- 
leading terms in the disjoining pressure. If the van der 
Waals forces direct the droplet towards the less wettable 
material, the droplet will stop as soon as the advancing 
contact line hits the step. Otherwise it will continue at 
a velocity rapidly decreasing with the distance from the 
step. Droplets which span the chemical step will move 
towards the more wettable substrate; however, the reced- 
ing contact line can be pinned by the step. 

This study demonstrates that taking into account the 
effect of long-ranged intermolecular forces is mandatory 
for accurately controlling and guiding the liquids in open 
nanofluidic systems. Recent experiments have shown 
that the arrangement of droplets on structured substrates 
can be explained by their crossing of chemical steps from 
the less wettable to the more wettable side Our 
study indicates that in general there can be free-energetic 
barriers to this process which would result in significantly 
altered patterns. Our analysis provides a microscopic 
approach to the pinning and depinning of three-phase 
contact lines at chemical surface heterogeneities which 
goes beyond the macroscopic picture of a sharp transi- 
tion between reg ions of different wettability on a sub- 
strate m, [2^ [so, [U, [111 or the phenomenological meso- 
scopic approach of introducing lateral variations of the 
jjax ameters entering into the effective interface potential 
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